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We present preliminary results of a study of the two-body _B-meson decays to a charmonium 
state (cc) and a K+ or if*" (892) meson using a sample of about 349 fb ^ of data collected with the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. Here cc indicates either the ric 
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state, reconstructed in the KsK^n^ and K'^K 7r° decay channels, or the he state, reconstructed 
in its decay to r;c7. We measure B{B° r]cK''°) = (6.1 ± O.Sstat ± l.lsyst) x 10"*, B{B+ ~> 
hcK+) X B{hc ric-y) < 5.2 x 10"^ and Z3(B° hcK"") x B{hc -qc-y) < 2.41 x 10"'', at the 90% 
C.L. 

Submitted to the 2007 Europhysics Conference on High Energy Physics, Manchester, England. 
PACS numbers: 13.25.Gv, 13.25.Hw 



The B decays to S-wave charmonium states, like J/tp 
and 77c, have been observed to occur with large branch- 
ing fractions {B) of the order 10~^ [if. Experimen- 
tal study of B decays to singlet states of charmonium, 
such as rjc and h^, is more complicated than the B de- 
cays to triplet states, such as J/V', fpC^S) or xa, be- 
cause one cannot exploit the cleaner signature of final 
states including a lepton pair. In this document, we re- 
port measurements of the branching fraction for the fol- 
lowing decay modes: B° ricK*° , B° hcK*° and 
B+ hcK+ [2]. We also reconstruct the B+ r]cK+ 
decay to be used as a "control sample". The branch- 
ing fraction of B'^ — s- rj^K*^ is currently known with a 
40% uncertainty, (1.2 ± 0.5) x 10"^ [||, while B decays 
to he have never been observed. The Belle collabora- 
tion studied the decay B'^ h^K^ with he — > 77c7 and 
reported B{B+ r]c-fK+) < 3.8 x 10"^ at the 90% 
C.L. for an invariant mass of the rjcj pair in the range 
[3.47,3.57] GeV/c^ No other B+ or B° decay modes 
with he have been studied yet. The he meson has re- 
cently been discovered by the CLEO collaboration as 
a narrow peak in Tp(2S) —^ rjc'j'n''^ decays at a mass of 
3524.4 ±0.7 MeV/c2 [|, and this observation was con- 
firmed by the E835 collaboration 0. 

In the simplest approximation, B decays to a charmo- 
nium state and a K or K* meson arise from the quark- 
level process b ccs . The colorless current C7^(l — 75)0, 
which can create the S-wave states like 77c and J/'(/', can 
also create the P-wave state Xci- It cannot, however, 
create the O^"*", 2++ and 1"' states XcO, Xc2 and he- 
Therefore B decays to any of these three states have to 
be ascribed to more complex mechanisms, such as the 
interaction of two color-octet currents In this sce- 
nario, B decays to XcO, Xci or h^ are expected to occur 
as abundantly as those to Xci- -B decays to Xc\^'^*^ have 
branching fractions between 3.2 x lO"'^ and 4.9 x lO"'^ [if. 
The B^ — > XcoK~^ decay has indeed been observed with 
a branching fraction of (lAtoil) x 10"" [H- However B 
decays to XciK^*^ and hcK^*^ have not yet been observed 
and upper limits on their branching fractions slightly ex- 
ceed 10"^ ;i]. 

In this analysis we reconstruct the rjc in the A'g(— > 
7r+7r")_ft'^7r^ and K~^K~n^ decay modes, the he in its 
decay to r/c7, and the in the mode K*° K~^it~. 
The K'^K^tt^ and if+_ft'"7r° final states are manifes- 
tations of the same decay mode, KKtt: they are cho- 
sen because they are among the easiest ije decay modes 



to reconstruct and have a rather large branching frac- 
tion, B{rie KKtt) = (7.0 ± 1.2)% Q. The rjej decay 
of the he is chosen because it is expected to comprise 
about half of the total he decay width We mea- 
sure ratios of branching fractions with respect to that 
of B+ rieK+, (9.1 ± 1^) X 10"4, to cancel the 17% 
uncertainty on B{rie — > KKn). 

The data used in this analysis were collected with the 
BaBAR detector at the PEP-II e+e" storage rings, and 
correspond to about 349 fb~^ of integrated luminosity 
collected at the T(45') resonance, comprising 384 million 
BB pairs. The BaBAR detector is described elsewhere [8(. 
Momenta of charged particles are measured in a track- 
ing system consisting of a five-layer, double-sided sili- 
con vertex tracker and a 40-layer drift chamber (DCH), 
both in a 1.5-T solenoidal magnetic field. Identification 
of charged particles is provided by measurements of the 
energy loss in the tracking devices and by a ring-imaging 
Cherenkov detector. Photons are detected by a CsI(Tl) 
electromagnetic calorimeter (EMC). The BaBAR detector 
Monte Carlo (MC) simulation based on GEANT4 [i| is 
used to determine selection criteria and efficiencies. 

The event selection is optimized by maximizing the 
quantity Ns/\^Ns + Nb, where Ns {Nb) represents the 
number of signal (background) candidates surviving the 
selection. Ns is estimated on samples of simulated 
events, while Nb is extrapolated from regions far from 
the signals on data. Simulated signal events and data 
are normalized to each other using the available mea- 
surements for B ^ rjc decays and assuming B = \ x 10"*^ 
for B decays to he- 

We select events with BB pairs by requiring at least 
four charged tracks, the ratio of the second to the zeroth 
order Fox- Wolfram moment 10] to be less than 0.2, and 
the total energy of all the charged and neutral particles 
to be greater than 4.5 GeV. 

Charged pion and kaon candidates are reconstructed 
tracks having polar angles in the region 0.41 < 9 < 
2.54 rad, at least 12 hits in the DCH, a transverse momen- 
tum with respect to the beam direction larger than 100 
MeV/c, and a distance of closest approach to the beam 
spot smaller than 1.5 cm in the plane transverse to the 
beam axis and 10 cm along the beam axis. A K*^ can- 
didate is formed from a pair of oppositely charged kaon 
and pion candidates originating from a common vertex 
and having an invariant mass within 60 MeV/c^ of the 
nominal K*'^ mass 
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Photon candidates are energy deposits in the EMC in 
the polar angle range 0.32 < 9 < 2.44 rad that are not 
associated with charged tracks, have energy greater than 
100 MeV, and have a shower shape consistent with that 
of a photon. A 7r° ^ 77 candidate is formed from a pair 
of photon candidates with invariant mass in the range 
[115,150] MeV/c^ and energy greater than 400 MeV. The 
mass of such candidates is constrained to the nominal tt*^ 
mass [l| when subsequently computing kinematic quan- 
tities. 

A ifg —^ TT+TT^ candidate is formed from a pair of 
oppositely-charged tracks originating from a common 
vertex and having an invariant mass within 20 MeV/ of 
the K'^-meson mass. Its measured decay-length signifi- 
cance is required to exceed three standard deviations (a). 
The candidate is constrained to the nominal if" mass [1| . 

The B'^''^ r]cK^'*° candidates are formed by pairing 
a K*'^ or candidate, referred to as the primary kaon, 
and a K'^K^tt^ or K~^K~tt^ combination with invariant 
mass in the range [2.75, 3.35] GeV/ c^. The mass range in- 
cludes the J/ip resonance. The 5+'*^ rjc'^K^'*'-' candi- 
dates are formed by combining a K*^ or candidate, a 
photon with energy exceeding 250 MeV, and a K^K^n^ 
or K^K^tt'^ combination with invariant mass consistent 
with the rjc mass. We perform a vertex fit to the B can- 
didates and require the probability of the of the fit 
to exceed 0.002. We define two kinematic variables: the 
beam-energy substitued mass, tties = \/ -E-bcam ^ Pb ^^"^ 
where pb {Eb) is the reconstructed B 
momentum (energy) and -Ebcam the beam energy, in the 
e+e" center-of-mass (cm.) frame. B candidates are re- 
tained if they have tties greater than 5.2 GeV/c^ and Ai? 
within [-24,30], [-40,30], [-34,30], and [-40,30] MeV for the 
KlK^TT^K*^'+, K+K-TT°K*°^+, K°K'^Tr^jK*°'+, and 
K'^K'n'^'yK*^'^ combinations, respectively. We also re- 
quire the absolute value of the cosine of the polar angle 
of the B candidate momentum vector in the e~^e~ cm. 
frame to be smaller than 0.9. 

To suppress background, the K^tt^ , K^K^, K^K^ 
and K^TT^TT^ combinations with invariant masses within 
30 MeV/c^ of the D°, and D+ meson masses Q, 
respectively, are excluded to form B candidates. We also 
remove K^K~ combinations containing a primary kaon 
where the invariant mass of the combination is within 
30 MeV/c^ of the cf) meson mass [l| . 

In events where more than one B candidate survives 
the selection, the one with the smallest JAiS] is retained. 
In cases of B candidates composed by the same final state 
particles, thus having the same value of JAiS], we retain 
the one for which the primary kaon has the largest mo- 
mentum in the e+e^ cm. system. 

The samples surviving the selection include a sig- 
nal component, a combinatorial background component 
given by random combinations of tracks and neutral clus- 
ters both from BB and continuum events, and a compo- 
nent due to B decays with a similar final state as the sig- 



nal. Such "peaking backgrounds" exhibit the same distri- 
bution as the signal in toes and AS, but their KKTT{'y) 
invariant-mass distribution (tox) is different. The sig- 
nal content on data is therefore obtained by means of 
a maximum likelihood fit to mx for all candidates hav- 
ing Toes in the signal region [5.274, 5.284] GeV/c^, after 
subtracting the combinatorial background. The mx dis- 
tribution for the combinatorial background events is ob- 
tained by extrapolating into the TOes signal region the 
mx distribution measured in the toes sideband, defined 
by Toes < 5-26 GeV/c^. The correlation between mx 
and Toes is found to be negligible in the relevant regions. 
A binned fit is then performed on the TOEs-sideband- 
subtracted mx distribution. 

We perform an unbinned maximum likelihood fit to 
the Toes distribution as follows. The B component, ac- 
counting for the sum of signal and peaking background, 
is modelled by a Gaussian function whose width is taken 
from the simulation and whose mean is fixed to the B- 
meson mass jl(] . The toes distribution of the combinato- 
rial background is represented by an ARGUS threshold 
function [ll|. The total number of events and the AR- 
GUS parameters are left free in the fit. The spectrum 
for candidates in the toes sideband is normalized to the 
Toes signal window by using the integrals of the ARGUS 
component in the two regions (Fig.[T]). 

In the case of 5° rjcK*^ ^ the TOEs-sideband- 
subtracted mx distribution is fitted to the sum of an 
rjc signal represented by a non-relativistic Breit-Wigner 
convolved with a Gaussian resolution function, a J/ip 
component modelled by a Gaussian of the same width, 
and a background component accounted for by a first- 
order polynomial with free coefficients. The masses of 
the rjc and the J/iji are fixed to the world average [H. 
The width of the r]c is fixed to the value measured by 
BaB^, 34.3 MeV/c^ [13]. The mass resolution modelled 
by the width of the Gaussian functions is fixed to the 
value determined on MC events, 11 MeV/c^. The num- 
ber of signal and background events is left free in the 
fit. We measure 185 ± 22 rjc and 59 ± 12 J/1I) candidates 
(Fig. [2]). The value for the fit is 54 for 56 degrees 
of freedom {Ndof)- Repeating the fit with no rjc com- 
ponent we get a /NdoF of 127/57. A similar fit of 
the 5+ VcK^ control sample yields 670 ± 32 rjc and 
149 ± 16 J/V' candidates with a x^/NdoF of 86/56. All 
of the values returned by these fits are summarized in 
Table D 



TABLE I: Number of rjc (Nr,^) and J/ip (Nj/^) events ob- 
tained from the fits described in the text with statistical errors 
only. 





Mode 






X^/NdoF 






185 ± 22 


59± 12 


54/56 


B+ 


{KKtv)K+ 


670 ± 32 


149 ± 16 


86/56 
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FIG. 2: Fit result (solid blue line) superimposed on the mES" 
sideband-subtracted mx distribution (points with error bars) for 
S° {KK-k)K*°. The (red) dashed line is the result of the fit 
with no r}c component. 



ing these fits with no component, we get a /NdoF 
of 45/40 and 48/40. 



FIG. 1: Left: mES distribution for B° ^ {KKtt)K*° (top), 
B+ ^ ■qc^K+ (middle) and -+ ■qc^K*'^ (bottom) candi- 
dates; points with error bars are data, the solid line represents 
the result of the fit described in the text, the dotted line rep- 
resents the ARGUS component. Right: invariant mass distri- 
bution of the KKtv{j) system for B° ^ {KK-k)K*° (top), 
B+ ■qc-iK+ (middle) and B° ^ ■qr.^K*^ (bottom) de- 
cays; points with error bars are data in the m-ES signal region, 
the shaded area represents the background expected from the 
rriES sideband. In the top plot, r^c and J/V" peaks are visible; 
peaking background events are also present as signalled by 
the excess of the data points above the shaded area outside 
the two peaks. No appreciable B component, neither signal 
nor peaking background, is observed for the B^ tJc'jK'^ 
and ric'yK*° cases. 



In the case of B+ ric^K+ and B° VclK*°, 
the TTiEs-sideband-subtracted mx distribution is fitted 
to the sum of an he signal modelled by a Gaussian, and 
a background represented by a first-order polynomial. 
The mass of the he is fixed to the CLEO measurement, 
3.524 GeV/c2 Q. The Gaussian resolution is fixed to 
the value determined on MC events, 16 MeV/c^ [l3 |. 
In the fit, the number of signal and background events 
is left free. The fit is performed over the mx range 
[3.3,3.7] GeV/c^. It yields 11 ± 6 and 21 ± 8 h^ can- 
didates with a x^/NdoF of 41/39 and 42/39 for the B+ 
and 5° yields, respectively (Fig.[2]and Table [TH) . Repeat- 



TABLE II: Number of he (iVhJ obtained from the fits de- 
scribed in the text with statistical errors only. 



Mode 




xVNdoF 


B+ ^ iVc'y)K+ 


11 ± 6 


41/39 


5° ^ (r,.7)^*° 


21 ± 8 


42/39 



The stability of the fit results is verified for various 
configurations of the fitting conditions. We float the rjc 
mass and width, which are poorly known, and the he 
mass. The values for the signal yields and the floated 
parameters returned by these flts are consistent with the 
nominal configuration. We verify the goodness of tlie fit 
with the chosen model using a MC technique: we sim- 
ulate a number of experiments by randomly generating 
samples of events distributed in mx according to the 
models used in the fit. The number of events generated 
is equal to the number of events in the corresponding real 
data sample. The parameters of the distributions are set 
to their fixed or fitted values. The fit is repeated in the 
same conditions as on real data. The pulls for the num- 
ber of signal and background events are distributed as 
expected. The robustness of the fit is tested on simulated 
events by varying the number of signal and background 
events input, including the null result. The number of 
events returned by the fit is consistent with the inputs 
for all cases. As an additional cross-check, we verify that 
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FIG. 3: Fit result (solid blue line) superimposed on the m^g- 
sideband-subtracted mx distribution (points with error bars) for 
_B+ rjc^K'^ (top) and fj" rjc^K*'^ (bottom). No significant 
he signal is evident. The (red) dashed line is the result of the fit 
with no signal component. 



the observed number of J/?/; candidates in the data agrees 
with the expectations. 

We evaluate systematic uncertainties on the number of 
signal candidates by individually varying the parameters 
that are fixed in the fits by iblcr from their nominal val- 
ues. We also estimate the systematic uncertainties that 
arise from a different choice of binning, fit range, and 
background parameterization. The large natural width 
of the rye introduces the possibility of interference effects 
with non-resonant B decays with the same final state par- 
ticles that can modify the mx distribution with respect 
to the one used in the fit. The fit is repeated including 



an interference term between the rjc and the background 
in the fitting functions. The amplitude and phase of the 
interference term are left free in the fit. The variation 
of the rjc yield with respect to the nominal fit is taken 
as a conservative estimate of the systematic error due 
to neglecting interference effects. The total systematic 
uncertainty on the signal yield determination, summing 
in quadrature all the contributions, is 8.1, 4.3, 24.8 and 
18.1% for B+ r]cK+, B" -> r]cK*°, B+ hcK+ and 
B° hcK*^, respectively. Being evaluated as yield vari- 
ations on the data, most of these systematic uncertainties 
should improve with larger statistics. 



The selection efficiency for _B+ 



VcK^ 



IS 



6%. 



The ratios of the selection efficiencies with respect to 
B^ rjcK^, estimated by using simulated events, are 
0.64 ± 0.01, 0.51 ± 0.01 and 0.29 ± 0.02 for B° r]cK*°, 
B+ hcK+ and B° hcK*°, respectively. Most un- 
certainties on the efficiency cancel out in the ratios be- 
cause of the similar final states. The remaining uncer- 
tainties are mainly due to differences between real data 
and simulation in the photon reconstruction as estimated 
from photon control samples on data (1.8%), and the un- 



known polarization for i?" hcK*^ estimated as in [14 1 
(6%). 

Using the signal efficiency computed on MC events, 
the signal yield observed on data, and the number of BB 
pairs in the data sample, we derive B{B^ rjcK^) x 
B{ric KKtt) = (7.5 ± 0.4stat) x 10"^ This is in agree- 
ment with the world average (6.4 ±1.4) x 10~^ 

We calculate the ratios of the branching fractions with 
respect to B{B~^ rjcK^) using the ratios of signal 
yields and efhciencies with respect to i?"*" ^ ricK~^, 
Rr = r(r(4S') -> B+B-)/T{r{4S) B°B°) = 1.026 ± 
0.032 1], and B{K*" K+tt-) = 2/3, and summing 
the uncertainties in quadrature. Table lllll summarizes 
the systematic uncertainties on the measurements. We 
obtain: 



B{B° 7jcK*°) 
B{B+ ^ 7?, K+) 



= 0.67±0.09,tat±0.07syst, 



and the 90% C.L. upper limits 
B{B+ KK+) X B{hc 



B{B+ ^ ry, K+) 
B{B" ^ h,K*") X B{h, ^ 77,7) 
B{B+ ^ rjc K+) 



< 0.058, 

< 0.26. 



These are determined assuming that the measurements 
follow a Gaussian distribution around the central value 
with standard deviation given by the total statistical plus 
systematic uncertainty. 

Using B{B+ f]cK+) = (9.1 ± 1.3) x lO"", we derive 

S(SO ^ *0) = (6.1 ± 0.8stat ± 0.6,y,t ± 0.9br) X 10-^ 

where the last error is from the uncertainty on B{B^ — > 
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TABLE III: Summary of the relative contributions to the sys- 
tematic error on ii^.x- = 13{B° ri^K*°)/B{B+ -q^ K+), 
Rh^K = B{B+ h,K+)xB{hc rjc'y)/B{B+ -> r/^ K+) and 
Rh^K' = 6(5" ^ h,K*'') X B{h, ^ r>cl)/B{B+ ^ rjc K+). 





RvaK' 


a{R)/R (%) 

Rh^K 


Rh^K* 


Signal yield extraction 


4.3 


24.8 


18.1 


Signal efficiency 


1.4 


2.2 


6.7 


TjcK'^ yield extraction 


8.1 


8.1 


8.1 


Rr 


3.1 




3.1 


Total 


9.8 


26.2 


21.1 
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rjcK^), and the 90% C.L. upper limits 

B{B+ hcK+) X B{hc iid) < 5.2 x 10-^ 
B{B° -> hcK*°) X B(hc r^d) < 2.41 x lO""*. 

Finally, we calculate 

g(i?° h,K*^) X Bjh, ^ rjd) ^ 

at the 90% C.L. 

In summary, we obtain a measurement of B{B^ — )■ 
rjcK*^) in agreement with, and greatly improving upon, 
the previous world average. We obtain an upper limit for 
B{B~^ — > hcK^) X B{hc Vcl) in agreement with the 
result obtained by the previous Belle measurement, and 
a first upper limit on B{B^ hcK*^) x B{hc Vcl)- 
The results confirm the suppression of production in 
B decays with respect to the S-wave rjc and the P-wave 
Xci and Xco states. All results are preliminary. 
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